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ABSTRACT 


The Rankine assumptions were used as a basis for developing 
equations for calculating active and passive earthpressures within 
a slope extending to infinity. The analysis considers: cohesive 
and noncohesive soils, the angle of internal friction of the soil, 
seepage forces caused by a ground water table that is parallel to 
the ground surface, and the plane on which the stresses are to be 
found may be at any inclination. 


INTRODUCTION 


Beginning with Rankine more than a century ago, scientists and engineers alike have 
sought to improve techniques for calculating soil pressures beneath the earth's surface. 
This capability has been important not only for predicting soil pressures against 
structures such as retaining walls and foundations, but also for predicting the depths 
and pressures at which soil masses will fail when disturbed by either man or nature. 


This paper, based on Rankine's classic assumptions and the contributions of others, 
offers some additional refinements for calculating soil pressures. The equations have 
particular relevance to much of the steep, mountainous lands managed by the USDA Forest 
Service, particularly the Idaho Batholith, where Intermountain Station engineers and 
scientists are studying slope stability and erosion problems. 


In 1860, Rankine presented the original theory for calculating lateral earth pres- 
sures within a slope of infinite extent. Rankine's development was based on a conjugate 
stress relationship between vertical stresses and stresses on a vertical plane and consid- 
ered a dry cohesionless material with either a sloping or a horizontal surface. During 
the latter part of the nineteenth century, the Mohr circle of stresses was presented for 
the graphic representation of the state of stress. The analytical work of Rankine was 
subsequently adapted to a graphical method using the Mohr circle of stresses. 


In 1915, Bell extended Rankine's work to include cohesive soils. Bell's graphical 
procedure made use of Mohr's stress circle, but in an inconvenient manner. Martin (1961) 
derived an analytical expression that extends Rankine's work to include cohesive soils, 
but it includes neither seepage conditions nor a variable plan of investigation. 


This paper extends Martin's analytical development to include stresses on any 
plane in the slope, and includes seepage stresses in the analysis. It is assumed in 
the following development that seepage stresses are caused by a water table parallel 
to the ground surface. A chronological summary of the developments based on Rankine's 
assumptions is shown in table 1. 


The graphical methods presently being used to determine earth pressures are tedious 
and time consuming to use because, unlike the analytical expressions developed herein, 
they cannot be processed by computer. 


Table 1.--Chronologtcal summary of the development of infinite slope earth pressures 


Factors Considered : Approach 
: : Angle of : Unit : Angle of the : : : 
Investi-; Angle : internal : cohesion of : plane of in- : Seep- : Graph- : Analyt- 
gator : of slope: friction : the soil: vestigation : age 2 1Cal 2 call 


Rankine X X vertical X 
Bell X X X vertical X 

Martin X X X vertical X 
Hartsog X X X any X X 


STRESS RELATIONSHIPS ON A 
PLANE PARALLEL TO THE GROUND 
SURFACE WITHIN AN INFINITE SLOPE 


The term, infinite slope, is defined as a slope of constant inclination of unlim- 
ited extent. For practical purposes, a slope that has a length to depth ratio of 20:1 
or greater has been considered an infinite slope. 


The classical infinite slope analysis requires that soil properties be constant at 
any and all vertical sections along the slope. If this requirement is met, the slope 
can be studied by analyzing the stress conditions within an elemental volume of soil. 


Figure 1 shows the basic assumed geometry of the infinite slope. In figure l, Z 

is the depth to the plane of investigation, Z, is the depth to the water table, and 
(Z-Zy) is the vertical distance from the water table to the plane of investigation. 
The total unit weight of the soil material above the water table is denoted as y,, and 
Ysat 1s the saturated unit weight of the soil. The angle the ground surface makes with 
the horizontal is denoted as 6. It is assumed that the planes representing the ground 
surface and the water table are parallel. 


Figure 1.--Profile of an tnfintte slope. 


ELEMENTAL VOLUME OF SOIL 


Figure 2.--An elemental volume of sotl, AA'BB', within an infinite slope. 


The slope profile with an elemental volume of soil (AA'BB') is shown in figure 2. 
If a unit dimension is assumed normal to the plane of figure 2, and if AA' is assumed 
to be of unit length, the resulting plane (which shall be referred to as. plane A-A') 
will be of unit area. The force, W, acting on plane A-A', is equal to the weight of 
the vertical column of soil above plane A-A'. 


The equation for the vertical force, W, is found by adding the products of | 
respective unit weights and volumes within the vertical column of soil: 


W= YZ yCOS8 + Year (2-4,,) e058: (1) 
Equation 1 is valid for the general case; that is, plane A-A' is below the ground 
water table. When plane A-A' is above the ground water table, the quantity (2-2) is 
zero and ZT is replaced with Z. 
The vertical force, W, can be resolved into components normal (N) and tangential 
(T) to plane A-A', as shown in figure 3. The expressions for the normal force and the 


tangential force are, therefore: 


= 2 es 2 2 
N Vp 2 yCOS @ + Monge ZJcos 6, (2) 
T = ¥,2,Simocos 6 + Yoat (274,,) Sinecoss. (3) 


The normal stress, o,, shown in figure 4, is equal in magnitude to the force N, 
because plane A-A' is of unit area. It follows, from equation 2, that: 


= 2 7 26. 4 
Os V2 y COS 6 + Yeap Z cos 8 (4) 


Figure 3,--Force 
components acting 
on plane A-A'. 


Similarly, the tangential stress, t,, is equal in magnitude to the force T. From 
equation 3, the equation for t, can be written as: 


E. = Vz yStnocosé ar Y sat (2-2) Sin8cosé. (5S) 


Because the analysis will be on the basis of effective stresses, it is necessary to 
determine the neutral stress, u. The neutral stress is found by use of a flow net, as 
shown in figure 5. The pressure head along an equipotential line is equal to 

(Z-Z,,) cos*6. The neutral stress on plane A-A' is the pressure head multiplied by the 
unit weight of water CY) 


= - 2 
u y,, (2 Z J) cos Ole (6) 


Figure 4,-- 
Stresses acting B! 
on plane A-A'. 


Note: hydraulic gradient i = sin@ 


Figure 5.--Flow net within an infinite slope. 


The effective stress, Oa is found by subtracting the neutral stress from the 
total stress: : 


—— 2 es 26) 2 = 26, 7 
oO, Vp 2yoos 8 + Voor ZJ cos YZ Z J cos (7) 


Because the buoyant unit weight, Vp> is equal to (Veta - Ves equation 7 can be 
written as: 


es 2 = Z 8 
ee 42 COS 8 + Y, (2 Zz) cos 8. (8) 


Because water cannot resist shear stress at negligible velocities, the shear stress, 

Ta, 1s not affected by the neutral stress. Equations 5 and 8 are the basic relation- 
ships for the effective stresses at a given depth and on a given plane of investigation. 
With this information, the stress condition (or, ta) may be adapted to a Mohr diagram. 


Stresses cannot be added vectorially because they are not vector quantities, but 
because they are acting on the unit area of plane A-A', the normal and tangential 
stresses can be used interchangeably with normal and tangential forces. The subsequent 
development makes use of this resolving of stresses in order to define obliquity angles 
(the angle between the resultant and the normal effective stress). 


Figure 6.--Stress 
relattonshtps 
for plane A-A' 
above the water 
table. 


There are three different conditions for which the obliquity angle can be defined: 


Case I. Plane A-A' above the water table. When plane A-A' is above the water 
table the angle of obliquity is equal to the slope angle, 8, as shown in figure 6. 


Case II. Plane A-A' below the water table with the water table coinetding wtth 
the ground surface. The obliquity angle, 8, is defined in terms of effective stresses 
as shown in figure 7. The tangent of the obliquity angle is found by dividing the 
tangential stress by the effective normal stress: 


Y (Z-Z )sinécosé 
eames sat W ‘ (3) 


= 2 
(Z Z J cos ) 


Solving equation 9 for the angle of obliquity, 8, yields: 


F 
[(1+*) tang]. (10) 
b 


RESULTANT 


Figure 7.--Stress rela- 
ttonships for plane 
A-A' below the water 
table with the water 
table cotnetding wtth 
the ground surface. 


Case III. Plane A-A' below the water table with the water table below the ground 
surface. For this general case, shown in figure 8, the tangent of the obliquity angle, 


y, is the quotient of the tangential stress and the effective normal stress: 


‘n6 : ; 
Y,2,5in cos6 + Year 2 By eerecos 


tan wv = 
2 = 2 
Y,2 cos+é + Yp 2 Z )cos48 


Rewriting equation 11 and solving for the obliquity angle yields: 


_ =i Jil e2y* “eae 2e4) 
p = tan > | tani 
V2, + Y, (2-2) 


STRESSES CAUSED BY 
SOIL ABOVE THE WATER 
TABLE 


RESULTANT 


STRESSSS CAUSED BY 


TABLE 


(D) 


Figure 8.--Stress relattonshtps for plane A-A' below the water table 


with the water table below the ground surface. 


SOIL BELOW THE WATER: 


(11) 


(12) 


MOHR’S CIRCLE OF STRESSES 


The Mohr circle of stresses is used to find the state of stress on any plane per- 
pendicular to the plane of elemental volume of soil, AA'BB', shown in figure 2. Sign 
conventions used are those commonly used in soil mechanics; normal stresses in compres- 
sion are positive and counterclockwise shear stresses are positive. 


When the effective stresses, om and t,, are plotted on Mohr's coordinates, the 
obliquities of the various components are as shown in figure 9. Note that the stress 
condition on plane A-A' is a function of depth. The state of stress (Ga; T,) will lie 
on line DE if plane A-A' is above the water table, or on line EF if plane A-A' is below 
the water table. The state of stress (o,, ) shown in figure 9 can be adapted to a 


Mohr diagram, as shown in figure 10. 


(+) 


D STRESSES CAUSED BY SOIL ABOVE THE STRESSES CAUSED BY o 
WATER TABLE SOIL BELOW THE 
WATER TASLE 


Figure 9.--Stresses on plane A-A' shown on Mohr's coordinates. 


(5, nae b) pass 


(Ge yt e) pass 


fe) 


ACTIVE STRESS CIRCLE 


Qi 


PASSIVE STRESS CIRCLE 


N! 


Figure 10.--Mohr dtagran for active and passtve stress condittons. 


The strength characteristics of a soil material may be defined by Coulomb's 
equat. on: 


Ss = c + otand. (13) 
Where: s = the shear strength of the soil material; 
c = the unit cohesion; 
o = the normal effective stress; and 
> = the angle of internal friction. 
Equation 13 is shown graphically in figure 10, and this graphic illustration of equa- 
tion 13 is known as the Mohr rupture envelope. The leaves of the Mohr rupture envelope, 
MN and M'N', shown in figure 10, represent the strength of the soil based on Coulomb's 
equation. A state of stress above MN or below M'N' will cause failure of the slope, 
whereas a state of stress that lies within the envelope may exist in a stable slope. 
The leaves of the envelope define impending failure conditions. 
The main objective of this development is to find equations for two limiting 


values of stresses, commonly called the active and passive states of stress. When the 
soil is on verge of failure it will be in a state of plastic equilibrium. 
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ie 
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Figure 11.--Active state of stress. 


AL 


Figure 12.--Passive state of stress. 


In the active case the lateral pressure is relieved and elongation occurs parallel 
to the direction of pressure release, as shown in figure 1]. .The lateral pressure, P.,, 
shown in figure 11, is the minimum which will maintain stability. This state of stress 
is shown by the active stress circle in figure 10. 


In the passive case the lateral pressure is increased and compression occurs 
parallel to the direction of pressure increase, as shown in figure 12. The lateral 
pressure, P,,, shown in figure 12, is the maximum that will maintain stability. This 
state of stress is shown by the passive stress circle in figure 10. 


Figure 10 shows that three conditions must be met for a soil to be in an active or 
passive state of stress: (1) The stress circle must be tangent to the Mohr failure 
envelope, (2) the center of the stress circle must lie on the o axis, and (3) the stress 
circle must pass through the point of known stress (64574) 


A modification of the Mohr's circle of stresses will be used herein. This modifica- 
tion, the origin of planes, allows a better visualization of the stresses and the orien- 
tation of planes on which the stresses act than the conventional Mohr's circle method. 
The origin of planes, OP, is shown in figure 10 and is obtained as follows: 


1. Through the point of known stress condition, draw a line oriented identically 
to the plane on which the known stress acts. 


2. The intersection of this line and the circle of stresses locate the origin of 
planes (G>T) - 


3. The stress condition on any plane can be found by constructing a line oriented 
identically to the plane on which the stresses are to be found. 


4. The point where this line again intersects the circle of stresses yields the 
desired stress condition (Gote) 


The stresses for the active case (Gar hiee and the passive case (0 


) are 
Ge ~e7 pass 
shown for planes oriented 45° clockwise from the horizontal. : 
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ANALYTICAL DEVELOPMENT 


The problem is to develop expressions for the normal and tangential stresses for 
the active and passive states of stress, with the angle of the investigation, «, being 
one of the variables. Other variables in these expressions are: the known stress con- 
dition CART the slope inclination, 9; and the strength characteristics of the soil, 
c and >. There is no need to develop separate equations for (0,,1,) 4,4 and (osteo) pass 
because the development results in a quadratic equation that yields solutions for both 
the active and passive states of stress. 


_ The equation for a line parallel to the ground surface and passing through point 
(o >t), as shown in figure 13, is: 


t= 6 tan@é +b, (14) 


where b 1s the~ intercept. 


Substituting the known stress condition (o,,7,)5 given by equations 5 and 8, into equa- 
tion 14 yields: 


¥-2ySinbcosé + Y sat (2-4,) Sinbcos6 


(15) 
= Y,2,.cos*etane + Y, (Z-2,) cos*otane +b. 
Solving for b from equation 15: 
bes (Yat Vp) [(2Z-2,,) sinécos 6]. (16) 
Recalling that ee Oe equation 16 may be written as: 
br= ySin6(Z-Z )cos6. (17) 


Equation 17 is recognized as the equation for the seepage force per unit area, S: 


S = Yi = YSin6(Z-Z )cosé. (18) 
Where: 
i = siné (the hydraulic gradient as shown in figure 5); 
V = (Z-Z,)cos8 (the volume of the column of soil between the water table and plane 


A-A' shown in figure 2). 
The intercept, b, is therefore equal to the seepage force per unit area, S. Substi- 
tuting the expression for b, from equation 17, into equation 14, yields the equation 
for a line through point (ot); with a slope equal to tané: 


tT = otand + Yy(2-Z,)sinocosé = octane + S. (19) 


142 


T Q 


Figure 13.--Active and passtve stress circles and thetr relattonshtp to 
Mohr's rupture envelope. 


From equation 13 and figure 13, the equation for a leaf of the Mohr rupture 
envelope can be written as: 


t - otand = c = 0. (20) 


The normal form of equation 20 is the expression for all lines perpendicular to the 
leaf of the failure envelope: 


— + - ————_ = 0. (215) 
Yl+tan*$ Y1+tan7> Y1+tan“¢ 


Because the stress circle is tangent to the failure envelope, the radius of the stress 
circle, r, is perpendicular to the failure envelope at the point of tangency. Noting 
that t ) equals zero, the length of the radius is found by substituting the coordinates 
for the center of the circle of stresses (9 >To) into equation 21: 


a c 
ot ane 


Sy eee ee ee (22) 
Yl+tan*¢ VY1+tan2o 


The general equation for the circle of stresses passing through point Cre) aN 
therefore: 


re (o tang+c)? 
(o-0,)2 + 1% - ———_—_-—._ = 0. (23) 
l+tan2¢ 


The origin of planes, (Gist is defined by the intersection of the circle of stresses 
and the line parallel to the ground surface that passes through the point (Cat a)ie The 
general equation for this intersection is found by substituting the expression for Tt 
from equation 19 into equation 23: 
we 7 (a tang+c) 2 
(o-0,)* + (stan6+S)* - ——————— _ = 0. (24) 
l+tan*¢ 


Equation 24 can be rewritten to yield a quadratic ino: 


o2(1+tan26) + o(-20 ,+2tanés) 


(G tang+c)* 
oo 82-2 ee Ge (25) 
2 1+tan7¢ 


Equation 25 will lead to two solutions for co: one solution (dq) will correspond to 
the point of known stress CO xsta)s and the other solution (op) will correspond to the 
origin of planes (Tp) - Equation 25 can be solved by means of a "sum of roots'"! 
solution (Rosenbach and others 1958), with the following result: 


20 ,-2tan8S (26) 


1+tan26 


By substituting o, for 5, equation 24 can be rewritten as a quadratic in o 
which 6, is the only unknown: 


ne in 


Ae. « (ada 20 Zaps 
eae: oC 20, 2o tan o-2tandoc) 


+ o%(1+tan?$+tan7e+tan*¢tan76) (27) 


+ 26 ,taneS(1+tan“$) + S2(l+tan2o) - c% = 0. 
Solving equation 27 by means of the ''quadratic equation" and rearranging, yields: 
o5 = 5, (1+tan*9) + tanodc + [-o, (1+tan*9) -tangc]* 


- [52 (1+tan29+tan?e+tan?otan6) (28) 


1 
+ 20 tandS(1+tan*$) + $2(1+tan2$) - c?]!? 


The positive root defines the center of the passive stress circle and the negative root 
defines the center of the active stress circle. 
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Solving for 6, in equation 26, yields: 


(20, = 2eaneS 


(29) 
(1+tan26) 


Where: Gin is defined by equation 28; 
S is defined by equation 18; and 


o, is defined by equation 8. 


Because 6, has two roots, it follows that 6, will also have two roots, one for the 
active origin of planes and the other for the passive origin of planes, as shown in 
figure 13. An equation for 1 can be developed, in terms of the defined quantity, Gp, 
by substituting Oy into equation 19: 


Ty = Optane + y,(Z-Z,)sin@cos® = o,tand + S. (30) 
With the expressions, equations 29 and 30, defining the origin of planes, the state of 


stress, (Gg;7 ), on any plane of investigation, «, can now be obtained. The conventions 
and notation for the typical stress circle are shown in figure 14. 


(+) 


T 


Q 


Ql 


Figure 14.--Typtcal stress ctrcle. 


nS 


_ The general equation for the slope of the line, GH, that passes through point 
(op,Tp), as shown in figure 14, is: 


T=7 
b 

tans = ——.- (31) 
o-O, 


Therefore, the equation of the line that passes through point Capitals and has a slope 
of =, is: 


+ (s-o, )tane . (32) 


The state of stress on a plane that is oriented at an angle < can be obtained by 
developing an equation for the intersection of the stress circle'and line GH. The 
equation for the intersection is found by substituting the expression for t from equa- 
tion 32 into equation 23: 


_ _ (o tand+c)? 
(c=0.)7 + [t) + (s-o, ) tan=]* - —————— = 0. (33) 
l+tan29 


Equation 33 can be rewritten to yield a quadratic in o: 


o2(1+tan2<) + 5(-20 +21 tane-26, tan*«) 


b 


ee Pein Pa 
+ (ont, 2t) 5, tan *+op tan ) (34) 


(o tandtc) * 
See a eee set ie 
l+tan2¢ 


The two roots of equation 34 are determined by the "sum of roots" technique. One root 
is the normal stress at the origin of planes, Os while the other root is the normal 
stress on the plane of investigation, Oe 


25 _-2t, tan«+2G, tan2« 
Sy. 5-2, tane+26, = 


(35) 
€ l+tan2« 


b ’ 


Because all terms on the right-hand side of equation 35 have been previously defined, 
the only unknown, Gs is expressed in terms of defined parameters. An alternate form 
of equation 35 can be written by substituting the expression for as from equation 30 

into equation 35: 


20, + o,, (-1+tan?«-2tan«tané) -2taneS 


7 36 
a (36) 
: l+tan2« 
An equation for T, can be developed by substituting Sos into equation 32: 
Pak (6-0, ) tane. (37) 


Equations 36 and 37 are general equations for calculating the normal and tangential 
stresses on a plane at any inclination, when the soil within an infinite slope is at 
the active or passive state of stress. 
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EQUATION VERIFICATION 


Equation Compatibility 


To verify this analytical development, the equations presented will be compared to 
those developed by other investigators. Martin (1961) has shown that his solution for 
active and passive pressures on a vertical plane agrees with developments by various 
other investigators. Because this paper is an extension of Martin's work, demonstrat- 
ing that the equations presented herein agree with Martin's equations would also imply 
agreement with developments by other investigators. 


Compatibility is proved by showing that equation 35, for the normal stress G., 
reduces to Martin's equation for normal stress; therefore the restrictions (no seepage 
forces and stresses on a vertical plane) on Martin's development must be applied to 
equation 35. 


: ; ies : il ; 
Recalling the trigonometric identity, (1+tan*«<)=(——>..) equation 35 can be 
rewritten: 
2 


Gg = 6 c = ine fog 5. sinz< - o.. 
O25 20 cos 21) sin cos« + 20, sin oF (38) 


The normal stress on a vertical plane is found by setting « equal to 90° in equation 38: 


Oa 1G (39) 


Therefore, the normal stress on a vertical plane is equal to the normal stress at the 
origin of planes, oy in equation 29. op can be adapted to the no seepage force condi- 
tion by setting S equal to zero in equation 29: = 
20 

(@) - 


“pO ietanee as (40) 


Substituting the expression for Sg from equation 28 into equation 40, and recalling that 
Ge equals OF, for this special case: 


26. + 20 tan*p + 2ctand 
a a 
l+tan26 


- OG - 5 2 — 2 
Siege ap 


af 
(1+tan26) 2 


= 
x 


4[62(1+tan29+tan2o+tan?gtan26) - c?] 


(1+tan20) 2 


lL 


Equation 41 is exactly the same (except for notation) as Martin's equation for the 
normal stress on a vertical plane. Therefore, the general equations given in this 
analytical development agree with equations given by other investigators for their 
more restrictive conditions. 


Critical Depth 


The soil properties, slope of ground surface, depth of the water table, and depth 
to the plane of investigation will all affect the stress condition (Gas72) Active and 
passive stress circles can be drawn through point (0g,t,) if, and only if, this stress 
condition remains within the Mohr envelope. The Mohr diagram, with the stress condi- 
tions (65,73) superimposed, shows conditions where failure is impending. The critical 
depth can be defined as the depth at which the stress condition on plane A-A' is equal 
to the strength of the soil. When the stress condition (Gag ba) lies on the Mohr enve- 
lope, only one stress circle may be drawn. For this condition, the active and passive 
stress circles coincide, and failure is impending. 


The stress on plane A-A’', (G30); at the critical depth, in a cohesive material, 
is shown in figure 15. If: 


1. o>6>6;, no Critical depth exists; 
B>o>8, a critical depth exists; and 
Ss B2O>e, AcCrivical depth..ex1sts.; 


ies) 


The stress on plane A-A', (G54) at the critical depth, in a cohesionless mate- 
rial, is shown in figure 16. If: 


¢>8>8, no critical. depth. exists; 

o=B>6, no critical depth exists; 

¢=B=6, a failure condition exists at all depths; 
B>o>0, a critical depth exists; and 

B>6>o, a failure condition exists at all depths. 


ONBWNFR 


Figure 15.--Stress condition at the eritical depth when the unit coheston 
ts greater than zero. 
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Figure 16.--Stress conditton at the critical depth when the untt cohesion equals zero. 


The equations for critical depth are found by equating stress conditions to 
strength conditions. Equating the right-hand sides of equations 19 and 13, yields: 


otané + S = otand +c. (42) 


Substituting the expressions for on and S from equation 8 and 18 into equation 42 and 
solving for the critical depth, ae yields: 


Cre Y,2,,cos*6 (tané-tang) 


Zo 22 +——~——4 _________. (43) 
cr Ww Yp,C08 @(tanO-tand) + y sinécosé 
W 


For the special case of no ground water (Z-Z ) is zero. Solving equations 8 and 18 for 


Vth wee 
cr 


7 i al inte ete ee EN (44) 
cr y,cos*6 (tané-tang) 


which is the critical depth for the special case of no ground water. Note that in all 
the equations presented, a slope will be stable as long as the depth of the soil mate- 
rial is less than the critical depth. 


Figure 17 shows the stress condition at the critical depth for a typical slope. 


This figure is also helpful for visualizing the expressions presented in previous 
chapters for the stresses acting on plane A-A’. 


LS 


POINT DESCRIBING CRITICAL DEPTH 


Y,, (2-Z,,) sin@cos 


ee 
$ = \e 
AG 2 | 2 
¥, (Z-Z,,) cos Q Y,,(Z-Z,,) cos a 
2 


Y Z cosQsin@ 
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Figure 17.--Stress relattonshtps at critical depth. 
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NONDIMENSIONAL FORMULATION 


To reduce the number of variables, the equations presented in the analytical 
development were made dimensionless. Equation 36 for 6, was expanded and divided by 
c, the unit cohesion, to determine the basic dimensionless parameters: 


Voz 
LAY Se 
(e 
Y,, (Z-Z__) 
2 ee b W 
(e 


¥.(Z=-Z_.) Y 
at, <r Wawona us Gal 
Yb 


The development of the nondimensional expressions parallels the analytical develop- 
ment. It was necessary to define additional parameters for the equations that lead to 
the final nondimensional expressions for the normal and tangential stresses. Because 
each nondimensionalized equation has an analog in the analytical development, the 
equations are presented with few details. 


The nondimensional form of the effective normal stress on plane A-A', Cay is 
obtained from equation 8: 
oy 
N'! = — = Ycos?6 + Y'cos76. (45) 
a e 
The nondimensional form of the seepage force per unit area, S, is obtained from 
equation 18: 
S : YW F 
v =—= Xsinécosé = Y'{t{-——]sinOcosé. (46) 
Cc Vb 


The nondimensional form of the tangential stress on plane A-A', Ta is obtained from 
equation 5: 


1G 
TH = == = N"'tané + Xsin6coss. (47) 
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The nondimensional form of the center of the stress circle, Ces is obtained from 
equation 28: 


ro) 
noe _O _ ym 2 
No = Ni Ci+tan >) + tando 
+ [NiI2 (1+tan2$)? + 2N''(1+tan*9) tang + tan2o] 
-Ni'? (1+tan*o+tan*6+tan*$tan76) - 2N''vtan6(1+tan*¢) (48) 


1 
2 


-y* (l+tan2¢) +-1 


The nondimensional form of the normal stress at the origin of planes, Op» is obtained 
from equation 29: 


ep Oe 2 e _ Nw 49 
Ni cio = 2cos*@(N'-vtan6) - No. (49) 


The nondimensional form of the tangential stress at the origin of planes, Tad 
obtained from equation 30: 


fee, OPA " 
Ty = ; = Njtane Se ae (50) 


The nondimensional form of the normal stress on a plane at the depth of investigation 
oriented at angle -, Tos is obtained from equation 36: 
- 
Me cos*«[2N'4Ni'(-1+tan*«-2tan=tan6) - 2vtan«]. (51) 


The nondimensional form of the tangential stress on a plane at the depth of investiga- 
tion oriented at angle «, Tos is obtained from equation 37: 


T 


wef Lop Wont ee 52 
Le : (NY Nj) tan : (52) 


Figure 18 shows the notation on a nondimensional Mohr's diagram, which is 
analogous to figure 13. 
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Figure 18.--Nondimenstonal Mohr diagram. 


The equations for the nondimensional analogy of the critical depth were developed 
in the same manner as equations 42 through 44. The critical depth in terms of the 


nondimensional parameters is: 


1 - Xsin68cosé 
(Git Yo) = (53) 
— cos*6(tan6-tand) 


where: (Y + Y')¢, is the nondimensional form of the critical depth. For the special 
case of no ground water table, the nondimensional critical depth is: 


Y i pte Beet 2k eet (54) 
a cos*6(tané-tan¢) 
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TABULATED RESULTS 


Dimensional Form 


The equations presented in the analytical development were programed for solution 
on a digital computer. The slope inclination 6, the depth to the water table Zw» the 
depth of investigation Z, and the soil properties $, c, y,, and y,were assigned appro- 
priate values, and the active and passive stresses acting on a plane at inclination < 
were determined. The angle « was incremented from 0° to 360° to solve for the stress 
on selected planes. A listing of this computer program is shown in the Appendix, and 
typical results are shown in tables 2 and 3. 


Nondimensional Form 


The equations presented in the nondimensional formulation were programed for an 
IBM 1620, Model II, digital computer. The parameters 4, 6, Y, and Y', were assigned 
various levels, and the nondimensional values representing the active and passive 
stresses, acting on a plane at inclination <, were determined. The angle < was incre- 
mented from 0° to 360° in order to solve for the stress on selected planes. The 
parameters », 6, Y, and Y', were then systematically varied to provide a comprehensive 
tabulation of results. Typical results are presented in tables 4 and 5. 
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As an example, the resulting nondimensional tables may be used as follows; 


Lees ‘de="05 
Gna 20 
Ze 40) tt 
W 
Zama Ors) cet 


c = 320 lb/ft? 
= 41.60 1b/ft3 


Yb 
Yq = 180 lb/ft 3 
ie ase 


It is desired that the active stress for the above conditions be determined. The 
nondimensional parameters Y and Y' are calculated as follows: 


Yp(Z-Z,) (41.60) (6.3-4.0) 
YY! = —————— -=—@ —————_—--—————_ 2 0.5. 


c 320 


The upper portion of table 4 is applicable for the example values of ¢, @, Y, and Y'; 
therefore, values are obtained for parameters Ne and ss from the line corresponding to 
the «= value of 75°: 


N'' = -0.4902, 
Cc 


T'' = 0.5846. 
c 
oe 
Since NG es then To N= (320) (-0.4902) = -157.0 lb/ft2, which is the active 
i 
effective normal stress on a plane inclined at an angle « of 75°. Since ae = a then 


aS cr = (320) (0.5846) = 186.8 lb/ft2, which is the tangential stress on a plane 


inclined at an angle « of 75°. 
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CONCLUSIONS 


When the Rankine assumptions are to be used to determine earth pressures, the 
equations presented herein will yield the active and passive earth pressures, on a 
plane at any inclination, within an infinite slope having a water table that is parallel 
to the ground surface. 


The solutions presented do not include pressures from freezing or swelling of soil, 
hydrostatic water conditions, and so on. The solutions are documented for use when the 
Rankine assumptions are valid; however, the author does not advocate the use of 
Rankine's theory in cases where conjugate stresses do not exist. 
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APPENDIX 


HHKHHXKEXHXD TMENSTONAL COMPUTER PROGRAM## HHH H HX HHH 


STRESS JON ANY JPEANE *=TNPINIT TE SLOPE THEORY We HARTSOG 


CHKK KKK KKKHKERKEK EKER EH KK KEK AE KK KAKA HARKER AIK EK MAK HH KK HHH HH KERMA HH RH HK HK 


€ 
E 
€ 


iE 


READ ls» COHESs PHIDs GAMMA, THETAD 


COHES=COHESION PeSeFe 3 PHID=BRICTION ANGLE IN DEGREES 
GAMMA=UNIT WEIGHT OF SOIL AROVE WeTe IN PeCeFe 
THETAD=INCLINATION OF PLANE IN DEGREES 


l 


FORMAT (4F10e1) 
PRINT 29 COHES, PHID, GAMMA, THETAD 


ZOFORMAT (1LOXs1OHCOHESION =9F60194H PSF37Xs6HPHID =9F6el1 9511 DEGes 


17X93 7HGAMMA = 5F6els4H PCF /Xs8HTHETAD =9F60e195H DEGe) 


READ 3 »GAMMABs ZW 


GAMMAR=BUOYANT UNIT WETGHTs ZW=VERTICAL DEPTH TO WelTe 


3 


N 
33 


is 
34 


4 


72 


a= 


EORMAT (C2:F 1G. 1) 

READ 33%9N 
IS INCREMENT OF ALPHA IN DEGREES 

FORMAT (119) 

READ 345 Lo™M 
LS) INCREMENT OF DEPTH M IS MAXIMUM DEPTH 
FORMAT (2110) 

PRINT 49N»GAMMAB 9Z2W 

FORMAT (39X%93HN =s13915Xs8HGAMMAB =sF6els4H PCFsllXs4HZW =9F6els 
L4H Fle) 

PRINT 9- 

FORMAT (1HO 91 9X 9 46H x HHH KEKE HHH KK AK HK HK KKH K AHR KKH HK HI KH HH HH) 
PT=3e14159 

PHI=(PHID*¥PI)/180-6 

THETA= THETAD*¥PI/180. 

GAMMAW=62 .4 

A=GOSECTPHETLA) 

A2=A*A 

B=SINF( THETA) 

TANTSHB7 A 

TANT2=TANT#TANT 

H=SINF( PHI) /COSF( PHI) 

D=1e+H*H 

GA2=GAMMA#A2 

GA2D=GA2*D 

COHESH=COHES*H 

DO, Sah Ze 9/Mis ll 

Z=12 
LOFAL DEPTH TO PLANE UNDER STUDY 

LAC Z=ZW 2015; 20% BO 


20 21=2 


Z2=O6 
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GO TO 5 
30 721=7W 
C 21=DEPTH OF SOILS ABOVE W,T. 
22=2-ZN 
C 22=DEPTH OF SOILS BELOW W,T, 
5 S1G1=(21*GAMMA+Z2*GAMMARB) ¥A2 
C SIG1=NORMAL INTERGRANULAR STRESS ON PLANE INCLINED AT THETA 
S=B*GAMMAW*Z2*A 
TL=SIG1*TANT+S 
C Ti=SHEARING STRESS ON PLANE INCLINED AT THETA 
GY=S1G1*D+COHESH 
BOJ= ((-SI1G1*O-—COHESH) **2)-((SIG1*SIG1)*(D+TANT2+H*H¥TANT2) 
1+2e*SIG1L*TANT#S*#D4+S*S*D—-COHES*COHES) 
TE (BON) 109141992 41999 
1991 PRINT 199352 
19903 FORMAT (1H153H Z=5F5els68H*** FAILURE CONDITIONS EXIST AT THIS DEP 
1TH WITH ABOVE PROPERTIFS ***) 
1094 GO TO. 1005 
19092 BY=SQRTF(BOO) 
SI GOA=GO0-BO 
Cc. SIPGUA=CENTER -OF ACTIVE STRESS CYRCLE 
S1GOP=GO0+RO 
C SIGOP=CFNTER OF PASSIVE STRESS CIRCLE 
SIGZ2A=((26*SIGIA—2 eX TANT¥S)/¢1L0et+TANT2))-SIG1 
C SIG2A=NORMAL INTERGRANULAR STRESS AT ORIGIN OF PLANES FOR ACTIVE CASE 
S1G2P=((02.*S1GODP-2 e*TANT#S)/(1.+TANT2))-SIG1 
C SIG2Z2P=NORMAL INTERGRANULAR STRESS AT ORIGIN OF PLANFS FOR PASSIVE CASE 
T2A=SIG2A*TANT+HS 
C T2A=SHEAR AT OP. (ACTIVE CASE) 
T2P=SIG2P¥#TANT+#S 
C T2P=SHEAR AT OePe (PASSIVE CASE) 
C FOR ABOVE EQUATIONS REFER TO NOTES 
DFLA=ATANF((-SIGZA+SIGOA)/T2A) 
IF(DELA) 50559560 
59 TF (T2A) 59360560 
50 DDELA=1806F( 160./P1T)*DELAI 
C DDELA=ANGLE TANGENT MAKES WITH HORIZONTAL (ACTIVE) -NOTES 13-17 
GO TO 2901 
60 DDELA=(189e/PI)*DELA 
2901 DELP=ATANF((-SIG2P4+SIGOP)/T2P) 
IF(DELP) 1955209520C 
209 IF (T2P) 509449005409 
2090 TFIT2P) 390994094499 
490 DDELP=(180./PI)*DELP 
C DDELP=ANGLE TANGENT MAKES WITH HORIZONTAL (PASSIVE) -NOTES 13-17 
GO TO 80 
390 DDELP=180.+((180./PI)*DELP) 
GO TO g0 
199 TF (T2P) 50%,69%9699 
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500 DDELP=180.e+((180«e/PI)*DELP) 
GO TO 80 
600 DDELP=360e+((1806e/PI1)*DELP) 
80 PRINT 63 ZsDDELA 
6 FORMAT(1H153H Z=sF5elsl1Xs13HDELTA ACTIVE=oF 661) 
PRINT 7 
7 FORMAT (1HO/// 919X94HSIG1913X92HT1913Xs5HSIGOA911X s5HSIG2A912X%9 
13HT2A) 
PRIND 21s SIUGlsTlsSWGOAsSITGZAsTZA 
21 FORMAT (9X95Fl6e61l) 
PRINT 23 
DA FORMATO LHO 125 Xs SHALPHAs21X%s SHS IG3As12%o3H7 3A513X%s SHPSAy 
DO 24 IALPHA= No3605N 
AL PHAD=IALPHA 
ALPHAD=ANGLE OF PLANE FOR WHICH STRESS CONDITIONS ARE TO BE FOUND 
IF (ALPHAD-DDELA) 15915916 
16 ALPHAD=ALPHAD+1860. 
IF CALPHAD=369.) 15915292 
92 GO 10 10'5 
15 ALPHA =(ALPHAD*PI1)/180, 
TANA=SINF (ALPHA) /COSFL ALPHA) 
TANA2=TANA*TANA 
IF (ALPHAD-99.) 44940944 
44 IF (ALPHAD-2702) 48540548 
48 SIG3A=((2e*¥SIGOA-2 e® T2A¥®TANAH2 o¥®SIG2ZA¥TANA2)/( 1 0¢+TANA2) )-SIG2A 
SIG3A=NORMAL STRESS ON PLANE AT ANGLE ALPHA ~-{ACTIVE) 
T3A=T2A4+(SI1G3A-SIG2A)*TANA 
T3A=SHFAR STRESS ON PLANE AT ANGLE ALPHA -C(ACTIVE) 
GO TO 41 
BLOTS GS AH=Si1GZ A 
T3A=-T2A 
41 P3A=SQRTFISIG3A*SIG3A+T3ZA¥T3A) 
P3A=RESULTANT STRESS ON PLANE AT ANGLE ALPHA -( ACTIVE) 
24 PRINT 25 sALPHADsSIG3A9T3A9P3A 
25 FORMAT (1H0924X9F60199X%93F 1601) 
195 PRINT 669ZsDDELP 
66, FORMAT (His 3H Z=sFSe1,11¥sl4HDELTA PASSIVE=5F661) 
PRINT 79 
70 FORMAT (1HO///519Xs4HSIG1913X92HT1 913X95HSIGOP911X e5HSIG2P512X5 
13HT2P) 
PRINT 216 SiGlet bsSIGOP$S1G2P,1T2P 
PRINT 230 
230 FORMAT (1HOs25Xs5HALPHA 321X35HS1G3P512X 5 3HT 3P513X9 3HP3P) 
IF (SIGOP=SI1G2P) 81,82,82 
82 DO 78 IALPHA=N53699N 
AL PHAD=IALPHA 
IF (ALPHAD=DDELP) 11591159116 
116 ALPHAD=ALPHAD+1 80. 
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IF (ALPHAD-36%.) 11559115571 
11. GO TOS 
115 ALPHA =(ALPHAD*PI1)/189. 
TANA=SINF( ALPHA) /COSF( ALPHA) 
TANA2Z=TANA¥TANA 
IF (ALPHAD-90-0) 765374576 
76 IF (ALPHAD-270.) (3374373 
73 SIG3P=((2e¢*SIGOP—2 6*T2P¥TANA+2 6*SIG2P¥#TANA2)/(1Le+TANA2) )-SIG2P 
SIG3P=NORMAL STRESS ON PLANE AT ANGLE ALPHA -( PASSIVE) 
T3P=T2P+(SIG32P-SIG2P)*TANA 
T32P=SHFAR STRESS ON PLANE AT ANGLE ALPHA -(PASSIVF) 
GO TO 75 
74 STG3P=SIG2P 
T3P=-T2P 
75 P3P=SQRTF(SIG3P#SIG3P+T3P*T3P) 
P3P=RESULTANT STRESS ON PLANE AT ANGLE ALPHA -( PASSIVE) 
78 PRINT 25 sALPHADsSIG3PsT3P5P3P 
GO: <7-'0> 51 
81 DO 84 IALPHA=N53609N 
87 ALPHAD=IALPHA 
IF (ALPHAD+180.-DDELP) 849859385 
85 IF (ALPHAD-DDELP) 86,4,51551 
86 ALPHA=(ALPHAD¥PI)/180.6 
TANA=SINFC ALPHA) /COSF( ALPHA) 
TANA2Z=TANA¥TANA 
IF (ALPHAD-99.)91399391 
91° IF (ALPHAD-276.) 933599,93 
93 SIG3P=((26*SI1GOP=2 .¥T2P*T ANA42 « *SIGZ2P*TANAZ)/( 1 6+TANA2) J=S1G2P 
T3P=T2P+(SIG3P-SI1G2P)* TANA 
GO TO 94 
99 SIG3P=S1G2P 
T3P=-T2P 
94 P3P=SORTFISIG3P¥SITG3P+T3P#T3P) 
95 PRINT 253sALPHADSSIG3P3T3P3P3P 
84 CONTINUE 
51 CONT INUF 
1995 CALL EXIT 
END 
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Headquarters for the Intermountain Forest and 
Range Experiment Station are in Ogden, Utah. 
Field Research Work Units are maintained in: 


Boise, Idaho 

Bozeman, Montana (in cooperation with 
Montana State University) 

Logan, Utah (in cooperation with Utah 
State University) 

Missoula, Montana (in cooperation with 
University of Montana) 

Moscow, Idaho (in cooperation with the 
University of Idaho) 

Provo, Utah (in cooperation with Brigham 
Young University) 

Reno, Nevada (in cooperation with the 
University of Nevada) 
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